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ABSTRACT 


We calculate the cross-correlation of the Intensity fluctuations 
betv^een different frequencies and finite bandwidth effects on the 
intensity correlations based on the Markov approximation discussed in 
Paper I. Our results may be applied to quite general turbulence spectra 
for an extended turbulent medium. Calculations of the cross-correlation 
function and of finite bandwidth effects are explicitly carried out 
for both Gan sian and Kolmogorov turbulence spectra. The increases of 
the correlation scale of intensity fluctuations are different for these 
two spectra and the difference can be used to determine whether the 
interstellar turbulent medium has a Gaussian or a Kolmogorov spectrum. 



I . INTRODUCTION 


Frequency correlations of pulsar scintillations caused by the 
irregularities of interstellar medium have been observed (Scott and 
Collins 1968, Komesaroff et al. 1971, Salpeter 1969, Rickett 1969, 

Lang 19/1, and Sutton 1971) and the effect of finite receiver bandwidth 
on the observed scintillations has also been detected (Rickett 1969) . 

A theoretical interpretation of these phenomena requires a knowledge 
of the cross-correlation function of intensity fluctuations between 
different transverse coordinates and different frequencies, Pj(z, kj^, 
^ 2 ). (Lee 1974). 

Little (1968) and Lovelace (1970) calculated the effects of 
finite bandwidth on the observed intensity scintillations based on a 
"thin screen" model for a Gaussian irregularity spectrum. The "thin 
screen" approximation is not necessarily realistic in interstellar 
scintillations since the interstellar irregularities may not be confined 
to a small region between source and observer. Shishov (1975) calculated 
the frequency correlations for both the "thin screen" model and an 
extended medium. However, his results do not include finite bandwidth 
effects on the correlation function of intensity fluctuations and 
are limited to a turbulent medium with a single-scale power spectrum. 
Budden and Uscinski (1970, 1971) also calculated the effects of finite 
bandwidth on the scintillations of extended radio sources and found 
them to be practically negligible, in contrast to the conclusion 
of Little (1968) . 
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Tn Paper I of this series (Lee and Jbkipii 1975a), wo developed 
a theory of strong scintillations, which was applicable to a number 
of pulsar signals. Angular broadening and phase correlation were 
discussed in this first paper. In Paper II (Lee and Jokipii 1975b) 
we applied this theory to the temporal broadening of pulses. In Paper 
III (Lee and Jokipii 1975c) , a theory of intensity fluctuations for 
monochromatic waves was presented. In the present paper, we apply the 
techniques of Paper I to calculate the cross-correlations of intensity 
fluctuations between different wave-numbers and different transverse 
coordinates P^ (z, ki, ^ 2 ' ^ 2 )/ and finite bandwidth effects on 
the intensity correlatio. functions for an extended medium. The 
calculations are carried out for both Gaussian and power-law (Kolmogorov) 
spectra of the turbulent medium. We also point out the differences 
in finite bandwidth effects between a Gaussian spectrum and a power- 
law spectrum. This difference can be used to determine the nature 
of the interstellar irregularity spectrum. 
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II. GENERAL CONSIDERATIONS 


We consider the propagation for a wave E (r, t) whose 
Fourier component 

obeys the scalar wave equation 
2 2 

V 4',,<r) e (r* ♦ (r) = 0 (2) 

(»1 - (1) (D - 

Here (o)/2n) is the frequency of the Fourier component c is 

the speed of light, and e (r) is the refractive index in which the 

0 ) ' 

wave propagates. The refractive index c (r) is a random function and de- 

u - 

pends on both the position r and the wave frequency o). As an examole, we con- 
sider in this oaper the propagation of the high frequency waves with , the 

P 

plasma frequency of the medium, in the plasma medium. This applies to the oropa- 

gation of the radio waves in the ionosphere, in interplanetary space, or 

in the interstellar medium. If Ng is the electron density, then we have 

2 2 

e (r) e l“(a) /w ) (3) 

(U p 

and 

= 4ir N^ e^/m, (4) 

where e is the charge and m is the mass of an electron (Paper I). 

Now let N and e (r) vary randomly, and let < > denote an 

© 0 ) 

average over an ensemble of propagation volumes. Following Lee (1974), 
we define 


Ng (r) 

= <Ne 

(r)> + 


(r) 


= <c 

(r)> + 

5e 

(?) 

to 

(0 


to 


6(r) = 

- 4tt e 

2 

(r)/mc2 
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c (r) ■= 6c /<c > 3 (r)/k^ 

K - Ul (i) ' 

k = “ <c (1 - 

c “ c mti)^ ' 


r = (z, C) 
C = (x, y) 

and 


S = (K, k), 

where k is the wave-number and ^ is the transverse position. 

We also define 

= u {z, 5 , k) ( 5 ) 

for a plane wave propagating in the z direction. Then we have from 


Eq. (2) 


+ ^2 u(z. 5 » k) + 6 (r) u(z,c,k) = 0 , ( 6 ) 


where now 3 (r) is a wave-frequency Independent random variable with 


zero mean. 


As in Papers I, II, and III, we will consider here a situation 
in which each Fourier component be a plane wave propagates 

freely from the source along the +z direction until it strikes the 
turbulent medium at the plane z—0. The observers are situated at the 
plane z > 0. 

The intensity of the Fourier component with wave-number k and 
observed at position (z, g) is Kz,^, k) = |u(z, g, k)| 2 . The 
normalized cross-correlation function of intensity fluctuation between 
two Fourier components observed respectively at (z, and (z, t ) 
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can be written as 


Pj(z» ^1* ^1* ^2^ 


- '* ^1 ^ 

<IjXl2> 


- 1 


<U, U * U-> Ur 


•1 *^2 ^2 
<ui U 2 * > "■ 1» 


( 7 ) 


where = I (z, C^, k^) and = u(z, k^) for i = 1,2. 

In particular, if we set k^ = k^, becomes the spatial intensity 

corr/' iation function, which was discussed in Paper III. If we set c; = q , 

'i '2 

is the frequency correlation function of intensity fluctuations. 

The cross-correlation function in Eq. (7) is for detectors with 
zero bandwidth. We consider now the case of detectors with a finite 
bandwidth. Let fj^ (k) be the normalized intensity response of the 
bandpass, i.e. 




f (k) dk = 1 
K 

Then we have the observed intensity at (z, g) 

I(z, p =^u(z, C, k)j^ 

The correlation function of intensity fluctiiation measured with 


a bandwidth function f (k) is 

K 


(z, K) U 3 ^^ 3 Ug%f^<kpfK(kpdkidk 3 

u^*> fj^ (k~) dkj^ p 


( 8 ) 


where the parameter K in (z, K) denotes the characteristics 

of the bandwidth function f to be specified in Section IV. The 

K 


scintillation index (K) as a function of K is defined as 
m^j (K) = Qj (z, ^2 ~ . 


(9) 
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It is the purpose of this paper to calculate P^, and 
From Eqs. (7), (8) and (9), it is clear that in order to obtain 
these quantities, one has to calculate the fourth order correlation 
function !‘2 2 (z» Sj^, S2 S3, S^) = <u^^ Uj 03* U4*>- 

A transport equation for f2 2 with different wave-numbers has been given 
in Eq. (42) of Lee (1974). This transport equation is very complicated 
and will not be solved here. 

However, as was shown in Paper III, the fourth moment ^ 
the complex amplitude u with same wave numbers (monochromatic waves) 
is simply related to the second moments <n the strong scintillation 
region. It can also be shown that this relationship between the fourth 
and second moments of the field u is the same as for the case in which 
the complex amplitude u has a joint-normal distribution. Hence it is 
reasonable to assume that in the strong scintillation, the fourth moment 
P- of the complex field u with different wave-numbers is also related 

^ f 2 

to the first and second moments as if the complex u has a joint-normal 
distribution. Under this assumption, Eqs. (7) and (0) can be written 
respectively as 


and 



by noting that <u> and <Uj^U2> are negligibly small in the region of 
strong scintillations. Thus instead of solving the transport 
equation for the fourth moment f'2,2’ will determine Pj^ ^ (z, )t^, 

^2' ^ 2 ^ ~ ^'^1 which is the correlation between u^^ and u^* at 

different transverse coordinates and with different wave-numbers. 
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III. CROSS-CORREIATION FUNCTION OF INTENSITY FLUCTUATIONS 

A trarisport equation for j^(z, C 2 * *<2^ " ^^2*^ 

can be obtained from Eq. ( 6 ) under the “Quasi-optics" and "Markov random 
process" approximatioas (Lee, 1974 ) , and one Sas 


ar 


1,1 


3z 


' ^ > ’’l.l " 7 + i2 )AJO> 


It ^ 
'^2 


2A„ir,^-5,),r 
“^1 


( 10 ) 


2 


where 7^'^ = 7,^ , i = 1 , 2 , and 




^ X 


<0 (z, C^) 0 (z , 


( 11 ) 


The solution of Eq. (9) was discussed in Paper II in relation to 
calculations of the temporal broadening of pulses. Following Paper II, 
we define kj^ = k + -j", k 2 = k - and assume Ak is small and [Ak|<<k. 

We also normalize the intensity of incoming wave such that u(z = 0, k) = 1. 
Then for a statistically homogeneous medium, one can write (Paper II) 


Fi,i (z, c , k , k.) = r (z, C, Ak) F (z, Ak) 

- ii-ZZ U - R 

where ” ^2 


Ak 2 


F (z, Ak) = exp { - 
R 4k 


Ag(0) 2 } 


and Fjj (z, 5 , Ak) satisfies the following partial differential 


equation 


9F ^(z, C, Ak) ^ i^ y 2 


3z 


2 k' 


^ 2 k= 


IA„(0) - A (O) F„ = 0 


B 


0 ' 


D 


The initial condition for Fp is F^^ (z = 0) “ 1. F^^ refers to the 

"pure refraction" effect because it represents the result of 

different transit times of different ray patlis due to a varying index 

of refraction. F refers to the "diffraction effect" arising from 
D 

2 r, 2 

the diffraction term {^1- - _2 ) j. Eq. (10), or 7 ^ P- in Eq. (13) 

^1 K, Q Q 


( 12 ) 


( 12 a) 


(13) 
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In order to proceed further, we have to specify the two-point 
correlation of the fluctuating index Cj^(r) (or 6 (f)), from which 
the function A„ (C) appearing in Eq. (10) or Eq. (13) can be 

P - 

obtained through Eq. (11). Since the medium is statistically homogeneous, 
the correlation function <Gj^ (r) Cj^ (r')> depends only on |r - r’|. 

It is convenient here to wor)< with the spatial power spectrum Pj.(q) 
of the fluctuations, which is related to the correlation through 

Pg (q) « d^r < (X) (x + r ) > e ^ 5 • (14) 

We consider two types of power spectrum for (q) : a Gaussian spectrum 

and a power- law spectrum. For a Gaussian spectrum, Pj. (q) is written 
as 


"e <q) - 


(Gaussian) 


(15a) 


where L is the correlation scale of the fluccuations in refractive index. 
G 


For a power-law spectrum, P^(q) has the following form. 


2 

Pj. (q) = Bp e 

(1 + q^L 2) “/2 (Power-law) (15b) 

P 

where L is the correlation scale of index fluctuations and I is 
P P 

termed the cut-off (or inner) scale. Usually and equation (15b) 

for a “— 2 — is the Kolmogorov turbulence spectrum. Note that in each 
case B may be related to the local mean-square electron density fluctuations. 
We have 



- 128 TT j^-4 L 3 2^ 

G e G e 

Bp = 128 Ti ^-4 r(?)A(^ - | ) 

©2 _ - 
where r^ - is the classical electron radius and <^N 

IMC q 

fluctuating electron density. 


is the r.m.s 


(16a) 

(16b) 
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'or the Gaussian spectruffj of Eq. (15a) « we have from Cqs. (11) « (14) 


'« 


. »o k"! 

° 


-(- 


4Lq2 


‘ 41^ > (C«Lq). 


(17) 


Next we consider the power-law spectrum of Bq. (15b). For |(|>>S. , wc 

P 

neqlect the effect of the cutoff at q > and obtain froir* Bqs. (11), 

(14) and (15b) 


- Bpk4 U/Lp)>^ K^,(i;/L ) 

P Lp^ 22+U n r(p + 1) I 

whetx't pB^-l, 2p-f'|->0, and denotes a modified Bessel function 


of the second )cind. One can further show that for L c f 

,4 „ P 

t 1 - £12j;i) (^ ,a-2 j 

. 0 . 21^ > 


AbU) = 

^ 4if (a-2)L 2 

tr 


r(f) 


From Eq. (17) or(l8a), we may write 
Dg(C) « Ag(0) - Ag(0 

“ 0.x * 


(13a) 


(19) 


and thus Eq. (13) becomes 


ar. 


iAk ^2 


az ^ 2k^ 


5 


gp 

2k2 


= 0 


(20) 


For the Gaussian spectrum, v =■ 2 and 

,4 


6- 


Bg k‘ 


3^nL^ * 

For power-law spectrum with 2 < a < 4, 


B 


Bp k4 r ( 2 -S.) 


o 4t(a-2)2<^"2)r (“) , 


2 and 


(21a) 


(21b)' 
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Eq. (20) can be solved numerically (see Appendix, Paper II). 

The numerical results for Tjj, together with the values of in Eq. (12a), 
give us the values of which are then inserted into Bqs. (6*) and (9*) 

to obtain the cross-correlation function ^ 2 ' ^2^ 

intensity correlation of finite bandwidth, Q_ (z, C. , 

Before presenting the numerical results, let us consider the "pure 
refraction” effect, (z,Ak). This effect can be neglected if 

A0(o)s|«l . (22) 

For a Gaussian spectrum, the condition in Eq. (22) is easily satisfied 
in interstellar scintillations and the pure refraction effect can be 
neglected. For a KoJ-.- ;nrov spectrum (a = iL) in interstellar scintillations, 

Eq. (22) is not satisfied and the "pure refraction" effect is important. 
However, as discussed in Paper II, the "pure refraction" effect is caused 
by variation of the optical path along the line of sight and will not 
be observed at a position essentially stationary with respect to the scattering 
medium. In observations carried out over a period of hours or days, therefore, 
an Earth-bound observer could hope to detect only the diffraction effect Thus 

we will neglect the "pure refraction" effect in the following presentation. 

In any case, this effect could easily be incorporated into the calculation 
of Pj or because of its simple analytic form in Eq. (12a) . 

The numerical results for P^ (z, Cj^, )cj^, plotted in Figures 

(1), (2) and (3). Figure (1) shows the frequency correlation of intensity 
fluctuations, 

Pj(Af) « Pj (z, k^, ^2 = Si* ^2 " ~ — ) (23) 

c 
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as a function of thci normalized frequency separation F ■ (Af/f ) hotween 
two detectors with frequency f^^ and *" respectively. Curve (1) 

in figure (1) is for a turbulent medium with a Kolmogr-rov spectrum and 
curve (2) for a Gaussian spectrum. f|^ is the "decorrelation frequency >" 
which is defined as the frequency difference of two observing channels 
beyond which the intensities measured at these two channels are essentially 
uncorrelated. For an observer at z, we have the decorrelation frequency 

to - 2.C k 2'v + 2)/v M- (V + 2»/v ,2 


For a Gaussian spectrum v * 2, f^^ and 8^ is given by Eq. (21a) , while 

for a Kolmogorov spectrum, " ‘I’# given by Eq. (21b). 

Note that the decorrelation frequency f^ is related to the characteristic 

time tjj for pulse broadening (Paper II) and the characteristic scattering 

angle 0 (Paper I) respectively by 
c 

2ir fp tp “ 1 (25a) 

and 


f a ■■ . j 

D 

C 


(25b) 

for both the Gaussian and Kolmogorov spectra. 

Figure.. (2) and {it show, respectively for the Gaussian and Kolmogorov 
spectra, the cross-correlation function Pj(C, F) “ Pj (Zf 

as a function of Q - for various values o^ normalized frequency 

separation F ■= Af/f^. The characteristic intensity correlation scale 

t shown in the figures is related to the decorrelation frequency f by 
c D 




k 0, 


(26) 


and Eq. (25b). can then be written in tennis of the cliaracteristic parameters 
in the interstellar medium by Eq. (24). We only note that z~0.5 
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li,, s,SN for the Gaussian spectrum and C_« 

for a Kolmogorov spectrum. As shown in Figures {?.) and (3), 

the spatial cross-correlation Pj becomes broader as the frequency separation 

F is increased The broadening of P^ is caused by the preferential removal 

of the h.’gh- frequency components of its Fourier spectrun« a result of the 

incoherence between waves of different frequencies generated by the 
2 

term Ak Tjj in Eq. (13). 
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IV. BANDWIDTH EFFECTS 


In order to calculate the intensity-correlation function Qj 

for a detector with finite bandwidth* it is necessary to specify 

the bandwidth function (k) . As an example* we consider a rectangular 

bandpass of frequency width Af^ {or Ak^ * 2v&f^/c) t 

k - <k< k + — * ^o. and zero otherwise. Let K = Af /f_ be the 

normalized bandwidth. The correlation function can be obtained by 

inserting the numerical results of F into Eq. (8 ) and integrating 

1*1 

over wavenumbers. Figures (4) - (6) show the results. The scintillation 

index m as a function of the parameter K is plotted in Figure (4). 

Z 

Cu (1) is for a Kolmogorov spectrum and curve (2) for Gaussian spectrum. 

i’he scintillation index m is the unity for small bandwidchs and falls to a 

Z 

value of 0.5 for K = 5. 

Figures (5) and (6) show the renormalized intensity correlation function 
Q(j{i;*K) = Qj (Cf K)/Qj(<;,K = 0) as a function of C = j - ^ 2 ! 

bandwidth K for the Gaussian and the Kolmogorov 

spectra respectively. Note that the intensity correlation scale 

increases as the bandwidth K increases. {See also Little 1968) . As in 

the case of the broadening of P^* this effect is due to the preferential 

removal of high-frequency components arising from diffraction-generated 

incoherence. 

Let us define ? (K) be the half-width of the function 0 (?;, K). 

h N 

Figure (7) shows the percentage increase of ? (K) relative to ^ (0) as 

h h 

a function of the bandwidth K. As shown in the Figure, the percentage increase 

of C (K) for a Kolmogorov spectruin is almost twice as large as that 
h 

for a Gaussian spectrum. At K = 10* the increase of for a Kolmogorov 
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spectrum is 'v 60% while for a Gaussian spectrum, u. e increase is 
30%. The interpretation for this difference can be found from 
the power spectra of index fluctuations shown in Figure (6). 

Curves (1) and (2) show respectively the Kolmogorov and the 
Gaussian spectra. The observed wave intensity fluctuations are 
- sused mostly by the power inside the rectangle abed in Figure (8), 
where the spatial frequency q is close to the inverse 

of the Fresnel scale. The components of higher spatial frequency 
in the intensity power spectrum, which are removed preferently 
because of the incoherence between different frequencies for a 
detector with finite bandwidth, arise largely from that portion of 
the power spectrum contained in the area cefg in Figure (8). Inside 
the shaded area, there is more power to be removed from a Kolmogorov 
spectrum than from a Gaussian spectrum. Hence the higher frequency 
power in intensity fluctuations is reduced more severely and the 
intensity correlation function tiecomes broader for a Kolmogorov 
spectrum than for a Gaussian spectrum. Thus the difference between 
a power-law spectrum and a Gaussian spectrum are amplified here by 
the finite bandwidth of a detector. The broadening of the half-width 
of Qj can be used to determine whether the electron density fluctuations 
in the interstellar medium have a Gaussian or a power-law spectrum. 

Finally we remark that an analysis of the increase of the 
intensity correlation scale has not yet been carried out for 
interstellar scintillations. The published data on interstellar 
scintillations are not sufficiently accurate to determine whether the 
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interstellar turbulent medium has a Gaussian or a Kolmogorov spectrum 
(Lee and Jokipii 1975 d ) , because most observed quantities are not 
very sensitive to the for of the power spectrum. However, the 
increase of the correlation scale due to finite bandwidth of a 
detector is relatively sensitive to the form of the interstellar 
power spectrum. Therefore, we suggest here measurement of the increase 
of intensity-correlation scale as a method to determine whether 
the power spectrum in the interstellar mediisn is Gaussian or 
Kolmogorov. 


15 



FIGURE CAPTIONS 


Figure (1) 


Figure (2) 


Figure (2) 


Figure (4) 


Figure (5) 


Figure (6) 


This figure shows the frequency correlation of 
intensity fluctuations (F) as a function of 
the normalized frequency separation F = fif/f^. 
fp is the decorrelation frequency. Curve (1) 
is the for Kolmogorov spectrum and curve (2) 
for the Gaussian spectrum. 

The cross-correlation function Pj a function 

of q - ^ 2 ! for various values of the normalized 

frequency separation F"Af/fjj. is the characteristic 

correlation scale. This figure is for a 
Gaussian turbulent spectrum. 

As in Figure (2) for the medi'jm with a Kolmogorov 
spectrum. 

The scintillation index m_ (K) is plotted as a 

2 

function of the parameter K (normalized bandwidth) . 
Curves (1) and (2) are respectively for a Kolmogorov 
and a Gaussian spectrum. 

This figure shows the renormalized intensity correlation 
function (CrK) as a function of C = - C 2 I f°>^ 

various values of the normalized bandwidth K for 
a medium with a Gaussian spectrum. 

As in Figure (5) for the medium with a Kolmogorov 
spectrum. 
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Figure (7) 


This figure show.3 the percentage increases of 


the half-width C. (K) of intensity correlation 
n 

relative to Cu (K « 0) as a function of the 
n 

bandwidth K for both Gaussian and Kolmogorov 
spectra. 

Figure (6) This figure shows the power spectrum of the 

refractive index, (q) . Curve (1) is a 
Kolmogorov spectrum and curve (2) is a Gaussian 
spectrum. The intensity fluctuations of radio 
waves are caused mostly by the power inside 
the rectangle abed. The components of highest 
spatial frequency in the intensity power spectrum 
are mostly due to the power in the shaded area 
cefg. This figure is plotted in a log-log scaie. 
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Figure (7) 
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